Charge-memory polaron effect in molecular junctions 
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The charge-memory effect, bistability and switching between charged and neutral states of a 
molecular junction, as observed in recent STM experiments, is considered within a minimal polaron 
model. We show that in the case of strong electron- vibron interaction the rate of spontaneous quan- 
tum switching between charged and neutral states is exponentially suppressed at zero bias voltage 
but can be tuned through a wide range of finite switching timescales upon changing the bias. We 
further find that, while junctions with symmetric voltage drop give rise to random switching at finite 
bias, asymmetric junctions exhibit hysteretic behavior enabling controlled switching. Lifetimes and 
charge-voltage curves are calculated by the master equation method for weak coupling to the leads 
and at stronger coupling by the equation-of-motion method for nonequilibrium Green functions. 



Memory effects and switching at the molecular scale 
are in the focus of present experimental and theoretical 
studies within molecular electronics 0, 0, H, 0, H, H, 0] • 
Beside stochastic switching in single-molecule junctions 
[H, recent STM experiments 0, S| show multistability 
of neutral and charged states of single metallic atoms 
coupled to a metallic substrate through a thin insulating 
ionic film. The switching was performed by the applica- 
tion of a finite voltage to the STM tip and was explained 
by the large ionic polarizability of the film Q . 

The coupling of a charge to the displacement of ions 
in the film can be treated as an electron-vibron inter- 
action. If the energy of the unoccupied electron level 
without electron-vibron interaction is eo, the occupied 
(charged) state of the interacting system will have the en- 
ergy ei = eo — e p , where e p is so-called polaron shift (or 
recombination energy). Neutral and charged (polaron) 
states correspond to local minimums of the potential en- 
ergy surface and are metastable, if the electron-vibron 
interaction is strong enough. Applying an external volt- 
age, one can change the state of this bistable system, an 
effect that is accompanied by hysteretic charge-voltage 
and current- voltage curves. In this approximation it is 
not necessary to include Coulomb interaction explicitly, 
though one can additionally incorporate charging effects. 

It was suggested [1, [i| that bistability between charged 
and neutral states can be accounted for in a single-level 
model, when one electron level is coupled to one vibra- 
tion (Fig [p The same problem was also considered in 
Rcfs. [lfj, however with the conclusion that quan- 
tum switching between bistable states results in telegraph 
noise at finite voltage rather than in a memory effect. In 
this Letter we show that there is no contradiction among 
these two pictures, taking into account the time-scale of 
the switching process. Indeed, the switching time r be- 
tween the two states of interest should be compared with 
the characteristic time of the external voltage sweeping, 
t s ~ V (t) I \dV (t) I dt) . For r ^> t s , quantum switching 
can be neglected and hysteresis can be observed, while 



in the opposite limit, t -c t s , the averaging removes the 
hysteresis. We calculate the charge-voltage curves and 
describe the full crossover between two regimes. 

The Hamiltonian of the single-level polaron model is 



H = (eo + eipo)d d + ujQo)a + X (a' + a) d'd 
+ ^2 [( £lfc + eipi)c\ k c ik + (Vik4 k d + /i.e.) 



(1) 



where the first three terms describe the free electron 
state, the free vibron of frequency ojq (h = 1) and the 
electron-vibron interaction. The further terms are the 
Hamiltonian of the leads and the tunneling coupling 
(i = L,R is the lead index, k labels electron states). 
The electrical potential tpo plays an important role in 
transport at finite bias voltages V = ifL — V-R between 
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Figure 1: (Color online) (a) The energy diagram of the single- 
level electron-vibron model, coupled to left and right lead (or 
tip and substrate in the case of STM). (b) Franck-Condon ma- 
trix elements Mom, for weak (g = 0.1, squares), intermediate 
(g = 1, triangles), and strong (g = 10, circles) interaction. 
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the left and right electrical potentials, (po describes the 
shift of the molecular level by the bias voltage and can 
be written as tp = tp R + r](cp L - <p R ), 77 G [0, 1] [II]- 

The coupling to the leads is characterized by the level- 
width function r.;(e) = 27r^ fe | V^fc | 2 <5(e — e^), where the 
coupling Vik is assumed to be energy-independent (wide- 
band limit). The full level broadening is given by the 
sum r = Tl + Tr. 

Consider first the case of very weak coupling to the 
leads, T <§C uJo,e p . Using the polaron (Lang-Firsov) [H, 
[hi . fig canonical transformation, it is easy to show that 
the eigenstates of the isolated system (T = 0) are 



£(* t -*KV)»^E|o) 

\fm\ 



with the energies 

A 2 A 2 

E nm = ein + uj a m, e\ = e , e P = — • 

Ul U! 



(2) 



(3) 



When the system is weakly coupled to the leads, the po- 
laron representation, Eqs. (|2l3p . is a convenient starting 
point, n denotes the number of electrons, while the quan- 
tum number to characterizes vibronic eigenstates, which 
are superpositions of states with different number of bare 
vibrons. The qualitative picture of the sequential tunnel- 
ing through a polaronic state is given in Fig.flja). Here 
the potential energies of the neutral and charged states 
are sketched as a function of the vibronic coordinate x. 
When the external voltage is applied, the energy levels 
are shifted depending on the asymmetry parameter 7/. It 
should be noted that this type of the energy diagram is 
quite general for charge-controlled bistable systems. 

In the sequential tunneling regime the master equation 
for the probability p nm (t) to find the system in one of the 
polaron eigenstates ^ can be written as [3, 17, lj| 



dp r , 



dt 



Epnn' 
rn.Tr. 
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r n , n n 

1 m'mPrim 



I v [p}. (4) 



Here the first term describes the tunneling transition into 
the state \n,m) and the second term the transition out 
of the state \n,m). I v [p] is the vibron scattering inte- 
gral describing the relaxation of vibrons to equilibrium. 
The transition rates T™ m , are found from the tunneling 
Hamiltonian (the last term in Eq. (fl])). Taking into ac- 
count all possible single-electron tunneling processes, we 
obtain the incoming and outgoing tunneling rates at zero 
bias voltage as 



Tnim' — E/ r*(^ lm ~~ Eo m i) \M mm '\ 2 f®(Ei m — Eom'), 
i=L.R 

(5) 



i=L,R 



x (1 - f?(E lm , - E 0m )) . (6) 




Figure 2: Inverse life-time (iT) _1 of the neutral state (thin 
solid line) and the charged state (thick gray solid line) as 
a function of X/ujq at to = \ 2 /2loq; and the same at eo = 
0.9A 2 /^o (dashed lines), kT = O.lwo- 



Here / (e) is the equilibrium Fermi function, and 



= 



■ exp 



> is the 



Franck-Condon matrix clement. It is symmetric in to— to' 
and can be calculated analytically. For to < to' it reads 



^ (-ff)'y^!e-3/y"'-")/ 2 

M m<m ,=^ — — — ntn | — — — — , (7) 



1=0 



l\(m — + to' — m)! 



where g = (X/luq) 2 is the Huang-Rhys factor [201 ] . 

One characteristic feature of these matrix elements in 
transport is so-called Franck-Condon blockade 18, Hij ]: 



in the case of strong electron-vibron interaction the tun- 
neling with small changes in m is suppressed exponen- 
tially, as illustrated in Fig. \V[b) for the matrix clement 



M „ 



-9/23I 



1/2 



Hence only tunneling through high- 



energy states is possible. This is also suppressed at low 
bias voltage and low temperature. 

Finally, the average charge is (n)(t) = ^2 m Pim, and 
the average current (from the left or right lead) reads 

Ji=L,R.{t) = e 52mm' ipimm'POm' ~ ^imm'Pirn') ■ 

To proceed further, we calculate the characteristic life 
times of the neutral and charged ground states. The life 
time T nm of the state |n,m) is given by the sum of the 
rates of all possible processes which change this state, 
T nm = J2n'm< ^m'm- As an example, calculating the life 
time of the neutral state |0,0), with an energy higher 
than the charged ground state |1,0), we find 

T Q0 = X X — ^00) |-Mmo| fi(El m — E 00 ). 

(8) 

For energy-independent Ti (the wide-band limit) we ob- 
tain the simple analytical expression 

„m / \2 \ 

roo = r E e ~ 9 —^ £ o - - + u m . (9) 
to! V wq / 
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Figure 3: Inverse life-time (Yr) _1 as a function of normal- 
ized voltage eV/uo for the asymmetric junction (77 = 0) at 
X/ljq — 5 and eo = X 2 /2luo for the neutral state (thin solid 
line), the charged state (thick gray solid line) and the same 
for the symmetric junction (rj — 0.5, dashed lines). Inset: 
random switching between bistable states (dashed line) and 
single switching into the stable state (full line) after a sudden 
change of the voltage. 



The corresponding expression for the life time of the 
charged state is (assuming that the equilibrium electro- 
chemical potential in the leads is zero) 

T w = rVe-^/° -e + — + u m . (10) 

The dependence of the tunneling rates (|9I10[) on the 
scaled electron-vibron interaction constant A/wo is shown 
in Fig. [21 It is clearly seen that at large values of A 
the tunneling from the neutral state to the charged state 
and vice versa is exponentially suppressed in comparison 
with the bare tunneling rate T. Hence both states are 
(meta)stable at low temperatures and zero voltage. 

Based on the experimental parameters of Ref. [2J , the 
charged ground state is assumed to be below the equilib- 
rium Fermi energy of the leads, while the neutral ground 
state is above it. In the experiments Q the observed re- 
laxation energy e p ks 2.4 eV leads to the parameter X/loo 
of the order 5 to 10. Thus the system is in the blockade 
regime at zero voltage, see Fig. 

Next we consider the other important question, 
whether fast switching between the two states is possible. 
At finite voltage the switching rates are 

p —g n m 

T oo = E —T~ t r ^° ( £1 + w ° rn " C 1 " 7 ') ey ) 
L — ' ml 

m 

+ T R f{e 1 +cj m + V eV)], (11) 

p — g n m 

T io = E - J- t r ^° + ^ m + C 1 - i) eV ) 

rn 

+ T R f (- £l + Lo m - V eV)] . (12) 

The voltage dependence of the inverse life time (tF) _1 
is shown in Fig. [3] for a junction with the same tunnel- 
ing coupling, Tl = but asymmetric electrical field 




eV/co 

Figure 4: Population of the neutral state as a function of nor- 
malized voltage eV/uio in the asymmetrical junction (77 = 0) 
at \/uo = 5 and eo = \ 2 /2loo for fast voltage sweep (thin solid 
line), slower sweep (thick gray solid line), and in the adiabatic 
limit (dashed line). Inset: sketch of voltage time-dependence. 



(rj = 0), as well as for the completely symmetric junc- 
tion (77 = 0.5). The results in Fig. |3] imply that in both 
cases one can tune (rr) _1 upon sweeping the bias voltage 
and thereby control the timescales for switching between 
charged and neutral states. For the symmetric junction 
both switching rates, Tqq 1 and t^ 1 , (dashed lines) are si- 
multaneously nonzero at finite voltage (eV/ujQ > 40 for 
the parameters of Fig. [3]) leading to random switching 
(noise) sketched as dashed line in the inset. On the con- 
trary, for the asymmetric junction controlled switching 
into the neutral (black solid line) and charged (grey line) 
state can be achieved at large enough negative and posi- 
tive voltage, respectively. This qualitatively different be- 
haviour is a result of the distinct voltage asymmetry of 
the two inverse lifetimes which are never both finite. The 
further peculiar feature of the asymmetric case, namely 
that the switching rates of the neutral and charged states 
interchange their role as a function of bias, i.e., the neu- 
tral (charged) state is long-lived at negative (positive) 
bias, implies hysteretic behavior and a memory effect. 

To this end we consider what happens, if one sweeps 
the voltage with different velocity (Fig. 0J) for the asym- 
metic case 77 = 0. If the voltage is changed fast enough, 
i.e. faster than the life time of charged and neutral states 
(t ^> t s as discussed in the introduction), then both 
states can be obtained at zero voltage (hysteresis). In 
the opposite (adiabatic) limit the change is so slow that 
the system relaxes into the equilibrium state, and the 
population-voltage curve is single-valued. Note that this 
controlled switching is possible only for asymmetric junc- 
tions for the reason given above. 

We finally compare the results with those of a fur- 
ther important limiting case, namely that the level width 
is finite (and possible finite dissipation of vibrons is 
taken into account). Then the master equation ap- 
proach can no longer be used, and we apply alternatively 
the nonequilibrium Green function technique. Follow- 
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We showed that bistability arises if quantum transitions 
between neutral and charged states involved are sup- 
pressed, e.g. due to Franck-Condon blockade. Different 
regimes, characterized by random mutual transitions and 
by single switching events into a stable configuration are 
identified. In the latter case controlled switching of the 
molecule is achieved by applying finite voltage pulses. 

We acknowledge fruitful discussions with J. Repp. 
This work was funded by the Deutsche Forschungsgc- 
mcinschaft within the Priority Program SPP 1243 and 
Collaborative Research Center SFB 689 (D.A.R.). 



Figure 5: Average number of electrons at Ti = Tr = 5u>o 
as a function of normalized voltage eV/u>o for the asymmet- 
ric junction, r\ = (thin solid line), and for the symmetric 
junction, r\ = 0.5 (dashed line), for X/luq = 5 and e = A 2 /lj . 



ing Refs. [2]], the average number of electrons is 
determined by the lesser Green function G< (ti — tz) = 
i (dt (t 2 )d(ti)) as (n) = -if G<(e)^. The calculation 
of the Green function is a nontrivial task even in the 
single-level model. It is simplified in the important limit 
of low vibron frequencies, ujq <C T < e p , where the 
Born-Oppcnhcimer approximation holds true. We used 
the equation-of-motion approach in this case. In Fig. [5] 
the charge-voltage dependence, obtained in the simplest 
mean- field approximation [1, Q, is shown. The lesser 
function is represented as G < (e) = iA(e)f(e), with the 
spectral and the distribution function 

2r 

/(e) = T L fl(e-e<p L )+r R f R (e-e<p R ) _ 

The result is qualitatively the same as in the sequential 
tunneling case: For electrically asymmetric junctions two 
stable states exist at zero bias (memory effect), which 
can be switched by the voltage. The current shows simi- 
lar hysteretic behaviour as a function of voltage. For the 
symmetric junction hysteresis is observed only at finite 
voltage (nonequilibrium bistability 0]). Hence, asym- 
metric junctions are again preferable for a memory effect. 

Finally we note that in the case < T wc considered 
the stationary problem only, assuming that the switching 
rate between the two metastable states is small (com- 
pared e.g. to T) at large X/cJq. The calculation of the 
life times of metastable states within the Green function 
appraoch and of dynamical effects arising from the com- 
petition between voltage sweeping and switching times, 
such as in Fig. [4j remains as a problem for the future. 

To conclude, we considered a charge-memory effect and 
switching phenomena in single-molecule junctions taking 
into account dynamical effects such as the interplay be- 
tween timcscalcs of voltage sweeping and switching rates. 
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